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Abstract 

Dyonic classical solutions of Yang-Mills theory are considered and the 
complete set of fermionic zero modes of these solutions are studied. Rep- 
resenting the QCD vacuum as a gas of dyons, one obtains chiral symmetry 
breaking due to zero modes similarly to the case of instantonic vacuum. 

1 Introduction 

Dyonic solutions (DS) are known already for a long time [1-4] and were mostly 
used in the context of the QCD vacuum at nonzero temperature [5]. Using 
a (singular) gauge transformation one can present DS in a form where its 
magnetic and electric fields are time-independent and the former displays a 
clear magnetic monopole character [6]. 

One can also visualize DS in the 't Hooft ansatz [7,8] as a coherent infinite 
chain of instantons placed equidistantlly along the time axis. 

The case of a finite number of instantons was also considered [9] and 
it was shown that the magnetic monopole held is acting over a distance R 
approximately equal to the length of instanton chain L. 

The long distance held of DS can in principle bring about confinement, 
and simple estimates [9] show that a stochastic ensemble of dyons, e.g. non- 
interacting g produce an area law of Wilson loop, in agreement with 
the idea of the dual Meissner effect [10]. 

Recently dyons appeared in N = 2 SUSY Yang-Mills model [11] and 
it was shown that the strong coupling region corresponds in the model to 
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the weak coupling theory of dyons and monopoles. Condensation of those 
produces confinement again in line with the dual Meissner effect [10]. It 
is instructive to compare properties of DS and those of instantons. It is 
known that in the dilute gas of instantons confinement does not occur [12], 
nevertheless instantons may play significant role in the QCD vacuum. 

One important property of instantons is that it gives rise to fermionic 
zero modes (FZM) [13] and this brings in a lot of interesting physical con- 
sequences, e.g. the famous 't Hooft's determinantal interaction [14] yields 
proton instability, instantonic solution of the U(l) problem [15], and finally 
chiral symmetry breaking (CSB) in the instanton gas or liquid model [16- 
17]. It is the latter property which was studied analytically in [18-21] and 
recent numerical analysis of instantonic model vacuum [22] suggests that 
many realistic properties of qq correlators in the model are due to a reason- 
able description of CSB. 

A natural question arises whether DS can also ensure CSB. To this end 
we study in this paper fermionic zero modes of DS. For some cases some FZM 
are known explicitly, e.g. the Jackiw-Rebbi zero modes [23]. In our case an 
infinite set of FZM was found long ago [24]. 

We present in this paper a concise and economic formalism for DS and 
explicit form of its FZM. We also find a correspondence between one of the 
FZM and its Jackiw-Rebbi form. This we do in the first 3 chapters. 

In the rest of the paper we exploit the FZM of DS to study the question 
of CSB. To this end we consider a gas of dyons and apply the technique of 
instanton gas model to derive the effective propagation for quarks in the gas 
of dyons. We specifically study the formation of chiral (constituent) quark 
mass and chiral condensate in the dyonic gas and compare results with those 
of instanton gas model. 

The plan of the paper is as follows. In the next Section we adopt useful 
notations and economic formalism to rewrite gauge field in the 't Hooft ansatz 
and the corresponding Dirac equation for the quark zero modes. In Section 
3 we define exact dyonic solutions and study its properties. In Section 4 we 
classify fermionic zero modes on a dyon, and find correspondence with the 
Jackiw-Rebbi solutions [23]. In Section 5 we define the gas of dyons and 
in Section 6 we modify the formalism of the instanton gas model [18,21] to 
study CSB in the gas of dyons and antidyons. The last section is devoted to 
the discussion of results and perspectives of dyonic gas model. 
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2 Notation and formalism 



We begin our consideration by setting the notation. We are working in 
the Euclidean space-time, with indices running from to 3. Summing over 
repeated indices is implied. Let us define matrices and as follows: 

r„ = (/, -iff) (i) 

f f, = (I, iB) 

where a are the Pauli matrices. The rj and fj symbols are then given by 

where all indices take 4 values (0,1,2,3) and i] 0fll/ = fj 0fll/ = 6^ (our notation 
differs from that of 't Hooft). Given a real 4- vector we may define the 
matrices 

a = a^T^ (3) 



a — a^Y ^ — a + 

and indeed we have 



aa = a 2 I (4) 

with a 2 = a^a^. In what follows we will be considering both colour and spin 
indices and for that purpose we will write a c or a s depending on which space 
the matrices are acting on. 

Fermionic fields ip aa (x) carry both colour (a) and spin indices (a). The 
representation of the Clifford algebra that we are considering is 

^ = ( r° T o ) ' 75 = ( I -°i ) • ^ 

In what follows we will concentrate on chiral fields 75 ^/>± = i^±- The Dirac 
equation and operator for positive chirality fields ip + reduces to 

= , b = (d,- iA,)v, s (6) 

and for those of negative chirality ip_ is 

= ; D = (d, - iA,)V, s , (7) 
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where is the vector potential matrix and subscript s in implies that 
SU(2) indices in T refer to spin, whereas those in - to colour. 

In what follows we will concentrate ourselves on vector potentials of 't 
Hooft ansatz form [7,8] 

A; = -fjavM . (8) 

Defining A° = d^f we may introduce the matrix in color space 

A= l - ( 9 ) 

which differs from written before by a term proportional to the unit 
matrix. In this form 't Hooft ansatz can be written as 

A, = V , c (d c f){- 1 -) (10) 

where all matrices are acting on color space. In terms of this matrix A^ we 
may introduce modified Dirac operators 

V = (3, - iA,)T, s = D — l -{d s f) = e^De-y (II) 

v = (d, - iA,)T, s = b- l -{b s f) = e ybe-y 

An important role in what follows in played by the matrix M 

^ = ^Ev,c = ^Ev (ie (12) 

more explicitly M aa;f3b = \ ^(r \ s ) a p(J. \ c )ab] which acts on the 2 X 2 dimen- 
sional space of chiral spinors. This matrix is self adjoint, unitary and satisfies 
the following relations. 

M 2 = 1 
a s M = Mil 

a c M = Ma s (13) 

= Mil 
a \l = Ma, 
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From M we can construct the projectors 

p + = ki - m) = - £ r^rv = \ E f „.r , c (u) 

Z A 1 A 1 

P_ = i(l + M) 

which project onto a 1 and 3 dimensional space respectively. These spaces 
correspond to the total spin zero and total spin one spaces of the product 
of the original spin-isospin one half spaces. The dimensionality follows from 
TrM = 2. In the V + space, projected by P + , all vectors are proportional to 
the unique element of the b which is taken to be unitary. With this 

notation the modified Dirac operators for 't Hooft ansatz are given by 

V = 3 S - M(dJ) (15) 

T> = d s -2P+(d c f) (16) 

We close this section by giving some relations which will useful later in finding 
solutions to the Dirac equations. We have 

u + = r\ c u + (17) 

which can be proven by realizing that (u + ) aa oc (cr 2 ) a „. 



3 Euclidean dyons and antidyons. 

Using the 't Hooft ansatz [7,8] formula one can compute the field strength 
tensor F^. Separating it into self-dual and anti-self-dual parts 

= ri^J-^ + fji^u^ (18) 



we get 



(") 
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T lc W^W- x (19) 



(+) = — W^dTirdW 
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where 

F = w s e ~' (20) 

For self-dual solutions lo\ ^ vanishes. As seen from Eq.(19), if F\ and F 2 are 
two solutions of lo\ ^ = 0, then F\ + F 2 is again a solution provided they 
do not have overlapping singularities. Hence one can obtain solutions by 
combining the pure-gauge solutions: 

F = \ (21) 

f " 2 



(x — x ) 2 

We have, thus the two main series of solutions: 

N o 2 

F = l + Y- (22) 

N p! 

F = T [{x _ l x(i y ) 2 ( 23 ) 

Now we will specialize on dyonic solutions. These follow from the second 
of the above series (Eq. (23)) for Pi = 1 and x^ = (0, 27m) ( we intro- 
duce dimensionless coordinates dividing x^ and x^ by some arbitrary scale 
parameter). 

One obtains the form found by Harrington and Sheppard [2] 

„ 1 sinh r 

2r (cosn r — cos t) 

where t = x and r = \x\. The corresponding A^ a is time dependent. There 
exists however a general singular gauge transformation [9], such that 

~ (7 Ti 

A, = U + (A, + id,)U; U = exp(i-^-9), (25) 

where 

W 

tan0 = W o [ W r ]~\ (26) 

r 
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which in case of infinite multiinstanton chain (24) yields time-independent 
solutions (Rossi solutions 



1 r 

A m = -Uabnb{—. 1) (27) 

r smhr 

A 0a = -n a (r cothr - 1) (28) 
r 



The field of the dyon is selfdual, 



Eka = B ka = S ak — + n a n k (-— + 9 l ) (29) 



where g is 



9 = (30) 
smhr 

For an anti-dyon solution, which is the gauge transform of the multi- 
antiinstanton solution, one has 

1 r 

Aa = -Uabnb{—. 1) (31) 

r smhr 

A 0a = — n a (r cothr - 1) (32) 
r 

The magnetic monopole number of this solution Q is opposite to that of 
the dyon. 

The multidyon solutions are obtained from the 't Hooft ansatz (10)-(23) 
when one takes several parallel lines of instanton centers We shall not 
be using those explicitly below. 



4 Fermionic zero modes on a dyon. 

The set of fermionic zero modes is easily obtained from our formulas in 
section 2, specifically (6), (9), (11), (16). Here we present the result (see the 
Appendix for details of derivation). The set is one-parametrical, with real 
parameter f3 : < f3 < 1 

= F i/2fi( F -i F W)) u+ (33) 
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where F = Fhs, given by (24), and 

oo if3n2ir 

= T (34) 

This set of solutions can be connected to a set of zero modes found in 
[24]; for more details see the Appendix . 
Summing up in (34) one obtains 

ip\ e~ ttf3 sinh((l — |/3|)r) + cost ■ sinh(r\f3\) + isint sinh(r\f3\) 
Zr cosnr — cost 

In particular for f3 = | one has 

, n u sinh^ ■ cos ^ , 
r [cosnr — cos t) 

Insertion of (36) and (24) into (33) yields 

m ii , sinh^ ,i , 

5> = - ( rff) 5 £ '- (:i7 > 

where U is the inverse of the gauge transformation that yields the time- 
independent field strength (Eq. (25)), and given by 



/ 2 i ■ 1 r ■/ -> ->\ t r 

U = \ (sin— cosh — \- i(o~ r n)cos—sinh—) . (38) 

V coshr - cos V 2 2 V ; 2 2 J y ' 

Hence, the f3 = | zero-mode can be gauged transformed to a time-independent 
solution. All other solutions (for f3 ^ |) are in general time dependent. 

For the monopole classical solution it was shown by Jackiw and Rebbi, 
that there exists at least one normalizable (in 3d) zero fermionic mode x, 
which can be written as [23] 

x = -(g(r) + h(r)((f s n))u + . (39) 
r 

The only normalizable solution has the form [23] h = , 

f r dr' 1 

g(r) = const r exp[ / — (^( r ') -^( r '))] (40) 

Jo r' 2 
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where 

A(r) = 1 ; B(r) = r cothr - 1 . (41) 

sinhr 

Doing integrals in (39) one obtains 

^(1-e-n 1 / 2 _ r ^(tanh r -)^ 2 / N 

g(r) = constVr ^ + e _ r ^ 3/2 e - = const — — (42) 

which is seen to give rise to the gauge transform ol Eq. (37), demonstrating 
that the Jackiw-Rebbi zero mode coincides (up to a gauge transformation) 
with our solution ?/^ = 2). 



5 Dyonic gas model of the QCD vacuum 

We assume that the QCD vacuum at zero temperature can be realistically 
modelled by the dilute gas of dyons and antidyons. We start with an ansatz 
for gluonic fields: 

N 

A,(x) = J2A^(x) + B,(x) (43) 

i = l 

where B^x) comprises perturbative and other possible nondyonic fields with 
prescribed gauge transformation 

B,(x) -+ U + (x)(B,(x) + i d,)U(x) (44) 

Each of A«( x) is the field of dyon or antidyon which is made of solutions 
(27-28;31-32) by some shift i?« and 0(4) rotat ion of the time direction 
and global color rotation fiW 

AW(x) = AW(x,i2W,o;W,ftW) (45) 

where for each unit 4 vector one can define 

r = [( x - Rf - ((x, - i?,K) 2 ] 1/2 (46) 

t = { Xjl - (47) 
so that (omitting superscripts i everywhere) 

A^(x) = n {t)+ A ll (r,t)n {t) (48) 
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and is given by (8), (27), (28). One can see in (47) that the homogeneous 
gauge transformation for A^} is ensured by the property of 0, that — > 
QU } fi + = U + Q + } and all construction in (43) is gauge covariant. 

At this point one must specify in which gauge A^\x) in (43) are summed 
up. It is tempting to use the gauge ( we shall call it the Rossi gauge [6]), 
where A^ is time-independent as in (27-28), representing static monopole 
solutions of the Prasad-Sommerfield type [3]. However, in this case the sum 
(43) leads to an infinite action, whenever all possible orientations can be 
present in the sum. 

We have found only one gauge (modulo global rotations) where the action 
of the superposition (43) is finite. This is actually the singular gauge of the 
original 't Hooft ansatz (8), 

(W^r = I I" 11 " (49) 
2r coshr — cost 

yielding for a standard (not shifted and not rotated) solution the form 

A sinhr sint 
A} a ' = e atk n k ( cothr H -) + b ia - 



r coshr — cost coshr — cost 

.ivi A , sinhr 

4i = ~n a (- - cothr + 50 

r coshr — cost 



At large distances solutions (50) behave as 



A w ~ - F w ~ — (51) 



and the same is true for the sum (43). 

Hence the total action of a gas of dyons of finite length is finite. To 
ensure that dyonic gas could be a realistic model of the QCD vacuum one 
must investigate the following points: 

1) to check that the classical interaction between (anti) dyons is weak 
enough at large distances, so that the dilute gas approximation could be 
reasonably justified. 

2) to prove the existence of the thermodynamic limit for the dyonic en- 
semble (42), i.e. that the total action of the (big) volume V4 is proportional 
to the volume, when it increases. 
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One must also prove that the free energy of the dd calculated with quan- 
tum corrections has a minimum at a finite (and dilute) density. 

All these points are now under investigation [25] and we shall not go into 
details, because of lack of space, but rather assume that the stable dilute 
dd gas indeed exists. The main question for us in the present paper is then 
whether the CSB occurs in such a gas. This point will be studied below in 
section 6. 



6 Chiral symmetry breaking in the dyonic 
gas 

For the gas made of equal number of dyons N + and antidyons iV_, N + = 
iV_ = y in the big volume V4 we assume that the thermodynamic limit exists 
for the total action and other extensive quantities like the free energy, when 
N — > 00, V4 — > 00 and jf is fixed and finite. 

We shall use for the dyonic gas the formalism similar to that exploited 
for the instanton gas by Diakonov and Petrov [18]. 

The main problem to be solved in this section is: given fermionic zero 
modes on each of dyons and antidyons; find the full quark Green's function 
for the dyonic gas with the superposition ansatz (43). 

To this end we make the same interpolating approximation for the one- 
dyon quark Green's function as in [18], i.e. in the exact spectral repre- 
sentation of i = 1, ...N 

containing all modes n = 1,2,. ..00, we keep only zero modes and re- 
place the nonzero-mode contribution by the free Green's function, since they 
coincide at large n ~ yjp* . 

Thus with So = (-iD(B) — im) -1 one has 

SV(x,y) = S (x,y)+ £ "' )(a) "' )+(y) (53) 

, — im 

zero modes 

One can see that S^ diverges as m — > 0, but we shall show however that 
the total Green's function is finite for m — > if N + = iV_. 
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Using (53) and (43) one derives the total Green's lunction as in [18,21] 
to be 

s = So -£««(*) f-^-) . y { : )+ (y) (54) 

iik \im + V J ^ 

n,m 

where upper indices i,k run over all dyon numbers, 1 < z, k < iV, while 
lower indices n, m run over all set of zero modes of the given dyon with the 
numbers i,k. We have also defined 

V; k m = / u$ + (x)i(d - iB)u^(x)d*x (55) 



We keep here the field B to make the formalism gauge invariant; in estimates 
we systematically put equal to zero. Note that u+ and u m in (55) should 
have opposite chiralities, hence V tk refer to dyon-antidyon (dd) or opposite 
(dd) transitions, otherwise V* fc is zero. 

One can introduce graphs as in [18] to describe each term in (54) as 
a propagation amplitude from a dyon i to a dyon k through scattering on 
many intermediate (anti)dyons centers, with scattering amplitude of each 
center (dyon) being and transition amplitude from center j (excited to 
the s-th level) to center / (excited to the r-th level) being VjJ . 

The lower indices are not the only new element in (54-55) as compared 
to the instanton gas model [18-19]. The zero modes depend also on the 
Lorentz orientation of dyon, in addition to the color orientation fiW and 
position i?W of the dyon, see Eq. (45). 

v${x) = nW«W(s - i?W,cuW) (56) 

Our next task is to compute the matrix elements of ( . 1 > ) lk fixing initial 

_ v im + V ' nm 

and final states and averaging over all coordinates of intermediate dyons. To 
this end we introduce as in [18] the amplitudes D tk m and P lk m for even and 
odd number of transitions V respectively 

"T 1 r,ik I r>(i) n(k) r-,(i) r-.(k) (i) (k) \ V° '/ 



In the definition (57) it is assumed that amplitudes of returns to the initial 
and final center i ar k are not included in D tk , P tk and should be added sepa- 
rately (which makes Eq.(57) not an equality, but rather a symbolic equation). 
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This amplitude ol the return to the center j we denote as 

A mn = D% n (R( j \R( j \rt j \rt j \u;( j \u;W) (58) 

Since in D 33 integration over all intermediate coordinates [R^ k \Ct^ k \u^) is 
done, A mn does not depend on R^ 3 \^ 3 \u^ and is a constant matrix. 

Taking into account any number of returns to the same center j, brings 
about a matrix e mn , defined as: 

£mn = — (1 - imA)^ n (59) 
m 

With its help the equations, connecting P and D can be written as follows 

PL = -—Vni— l d*^dtoWdwW\v«e m ,D* lm (60) 

im im ZV^ J i 

D nm = -^F 4 j d 4 R«Un( 3 UJ 3 ^-V:i,e sls P 3 ^ (61) 

As a next step we separate out the dependence of P } D } V on lower indices 
and on fiW,fi( fc ). To this end we consider zero-mode solutions u^\x) in the 
form of (A. 11) and make Fourier transform 

v${p) = J u^{x)e ipx d A x = e vpo2vn u {t) (p) (62) 

It is important that vM\p) does not depend on n altogether. Therefore with 
the help of (56) one has 

v^^^KJ^r^M^J 3 )) = J ^^-^XUp), (63) 

v^Jp) = e-^^-^u+ip^n+^i-p^uipi) (64) 
where p l = 3l Wi p, and is 0(4) rotation transforming time unit vector into 

We introduce now "amputated" amplitudes d, /, w instead of D } P, V as 
follows 
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and similarly for f(p\p k ); according to (64)) one has w(p\p k ) = —p 
Insertion of these definitions into Eqs. (60-61) yields 

where we have introduced 

u(p) = e +tpo2 ™u(p)e ns e- tp ^ s u + (p) (68) 

n,s 

One can see in (67-68) that / and d do not depend on rotations in p\p k and 
the integration over du l there acts only on v{p 3 ), so that with the definition 

v(p) = J v(pj)d(jj J (69) 

one obtains 

iNpvp 



d(p) = fff 2 (70) 

and f(p) is expressed through d via (66). The definition (58) can be used now 
to obtain the selfconsistency relation, taking into account that at m — > 0, 
A ~ -\ and therefore one has 

1 

as a result of insertion of (65) and (70) into (58) multiplied with e mn , one 
has 

= 2vyv, r d*p M\p) 
U ° N J (27r) 4 M 2 (p)+p 2 1 ' 

where we have defined the average number of zero modes per dyon — n , 

y 1/4 

n ~ -^j— , b is the internal scale parameter of dyons, 
We also introduced the chiral mass M(p) 
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where we used the fact that v is averaged over all directions and should be 
proportional to the unit matrix in Lorentz and color space. 

Eq.(72) goes over into the corresponding consistency relation for instan- 
tons [18] when n = 1 and matrix e becomes a number, while u(p) is the 
Fourier transform of the 'tHooft's zero mode [13]. 

The solution d(p) (70) assumes the knowledge of the matrix e nS} while 
the consistency relation (72) imposes only one condition. Therefore the 
strategy of solution is as follows. From (65) one finds A mn = D^ n = 
$ (27r) 4 e~ 27rip °( m ~") u + (p)d(p)u(p) through d(p). It clearly depends only on 
the modulus \m — n\. Then inverting the matrix A mn one finds e mn from 
(71). Finally from (68)-(69) one finds v(p) and inserts it into (70), defining 
d(p) The cycle is thus completed, and should be repeated till the convergence 
is achieved. 

One can also study another basis of zero modes, namely that of (34). 
In this case dependence on f3 can be also extracted, indeed 

u f(p) =^ e *u« = J2e 27Tnt{p °- p) u {t \p) = 

n n 

= ^6(p-p -k)uV(p)=6 V3 , po] uV( P ) (74) 

k 

where we have introduced notation £[/3 Po ], implying that f3 is in the interval 
[0,1] and 6 - function should be moderated first, introducing finite number 
of centers N in the dyon (52n=-N /2) anc ^ considering limit N — > oo at the 
end. 

In this way one obtains the same equations (65-67) for /, d if the new 
definitions are used, e.g. 

D%. = J P [27r)4 ^ o] u + (^)0+'%V)O fc u(/)V P o] (75) 
and where in (66)- (67) now v is defined as 

v{p) v{p) = J dLou{p)e [po ^ o] u + {p) (76) 

and 

£[po,Po] = J d P J d P'ti\P,P ] £ wti\P',Po] ( 77 ) 
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From (75) one deduces that A^/ = 5ppiA(f3) and hence also eppi is diagonal 
due to the relation 

i 

£ PP' A fJifl" = ~i 8 l3l3" ( 78 ) 



and is equal to 

e pp , = V £ (/3) = V-^ 1 ^) ( 79 ) 



i 

m 

with 



A(/3) = / ^ s wpo]y + (p)d( P Mp) (so) 



The system (70), (72), (76), (77) and (79-80) is now complete. 

We now proceed to write down the quark propagator (54) in terms of 
functions d, f and finally in terms of the chiral mass M(p) (73). 

Following the same procedure as in [18], one can rewrite (54) as 



ns , 



p 2 2V4 \im V 1 — zmA, 
X [ dn {t) dJ t) (uli ) (p } J t) )u { ; )+ (p } J t) )+ dyon <-> antidyon)- 



+ (81) 

In the appendix we write explicitly the structure of the dyon zero-modes 
density matrix u[ ! 'uW + . Submitting in (81) expression (65-67) and (70), (72) 
we finally obtain S(p) in the form 

p + iM(p) 

s{p) = T+F" (82) 

This form justifies the meaning of M(p) as a chiral mass, i.e. an effective 
mass of quark due to CSB. It coincides with the form of S(p) for the instanton 
gas [18], however the explicit expression for M(p) (73) differs. 
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The most remarkable leature ol (82) is the disappearance ol the massless 
pole from So(p). One should have in mind of course that the form (82) is 
gauge-noninvariant and obtained neglecting confinement. If one takes into 
account these effects, as in [21], the pole structure in (82) is supplemented 
by the area law due to the string between the given quark and an antiquark 
and the pole is never present in physical amplitudes. 

From (82) one can easily compute the chiral condensate: 

< qq >Mtnk.= -i < qq >Eud = i < trS(x,x) > = 

It is nonzero thus confirming the phenomenon of CSB in the dilute dyonic 
gas. 

7 Conclusions and discussions 

We have described two different infinite sets of quark zero modes for an 
(anti) dyon. We suggest to consider these zero modes as a driving mechanism 
for chiral symmetry breaking (CSB) in a balanced dyon-antidyon (dd) gas, 
similarly to the CSB mechanism in the instanton gas [18]. The composition 
of the dd gas is by itself a nontrivial problem because of the long-distance 
field of dyons. We have found one specific gauge where the superposition 
ansatz of the dd vector potentials provides a finite action. 

The resulting dd gas has classical Coulombic - type interaction, so that 
dyons and antidyons are not confined. Then according to the 't Hooft's 
principle [26] one may expect Confinement of quarks (and gluons). The 
possible role of dyons in the Confinement mechanism was also emphasized in 
Ref. [28]. The numerical studies of the Wilson loop in the dd gas are now in 
progress [25]. But the main emphasis of the present paper is on the possibility 
of CSB in the dd gas. Indeed, if the CSB is present in the dd gas then one may 
hope to have a model of the QCD vacuum, which obeys both properties of 
confinement and CSB. At the same time, it gives a possibility of explanation 
of the simultaneous disappearance of confinement and chiral symmetry above 
the critical temperature. A similar proposal made by one of us [27] attributes 
Confinement to a different kind of multiinstanton configurations. Our result 
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might well be extendible to this model as well, since the main ingredients of 
our formulas are quite general. In section 6 we have calculated the one-quark 
propagator in the balanced dd gas and demonstrated the appearance of the 
chiral quark mass. 

The selfconsistency equation (72) has the form similar to the case of 
the instanton gas [18], and goes over into the latter when the number of zero 
modes on each center (dyon or instanton) tends to unity. Also the form of the 
quark propagator is similar and it is clear that the quark-antiquark Green's 
functions can be obtained by the same method, as used for the instanton gas 
in [18]. 

From (82)- (83) (and the corresponding qq propagators) one can see that 
for Nj > 1 the original symmetry Ul(Nj) Cg) Ur(Nj) of the QCD Lagrangian 
is broken down to SUv(Nf) <E> £/y(l), while both SUa(Nj) and Ua(1) are 
broken. Thus the topological charges of dyons, creating fermion zero modes, 
produce a universal mechanism for breaking both Ua(1) and SUa(Nj), found 
for instantons by 'tHooft [15] and Diakonov and Petrov [18]. We refer to this 
universal mechanism as the topological mechanism of CSB. 

Many of features of the problem under study are left unclarified. 

First, and the most important is the detailed numerical study of the dd gas 
with classical interaction and quantum corrections to yield the equilibrium 
density of gas and resulting numbers for string tension and chiral condensate. 

This question is now under study. 

Second, the effective action for large Nj and the bosonization procedure is 
not done in the paper. Therefore, the effective action for Goldstone-Nambu 
bosons is lacking above. This is planned for future publications. 

The most part of this work was done while the second author (Yu.S.) 
has been visiting the theoretical department of the Universidad Autonoma 
de Madrid. This author is grateful for kind hospitality to the members of 
theoretical department and all the staff. He also wants to thank P. van Baal 
and G. 't Hooft for useful dicussions. 

He was also partly supported by the Russian Fund for Fundamental Re- 
search, Grant 93-02-14937. A. G-A acknowledges finantial support by the 
CICYT grants AEN93-0693 and the EC network CHRX-CT93-0132. 
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Appendix 

In the appendix we give a short derivation of the solution of the Dirac 
equation in the background field of a multi-instanton solution of the 't Hooft 
ansatz form (Eq. (8)). Then we will specialize to the case of the dyon 
solution. 

Consider the field tensor constructed from a gauge potential of 't 
Hooft 's form. One can separate it into a self-dual and anti-self-dual parts 
lo\ + ^ and lo\ ^ as shown in Eq. (18). In terms of these parts one can compute 
the product of the modified Dirac operators T> as follows: 

VV = V,V,-2iLotT ls (A.l) 

VV = V jl V jl -%uj-Y ls (A.2) 

For purely self-dual fields u\ ^ = and one sees that there are no negative 
chirallity zero-modes. Thus, we may simply consider the equation 

V x (+) (x) = (B s - 2P + (d c f))x (+) (x) = (A.3) 

Now let us write the following ansatz for X^ + K x ) 

x (+)( x ) = (d s <f>)u + . (A.4) 

The ansatz follows from acting with d s on Eq.A.3 and subsequently with the 
inverse of the 4-dim Laplace operator A -1 . Substituting the form A.4 into 
A.3., we get 

(A<f>)u + = 2P + ((d c f)(d s <f>)u + (A.5) 

Now using the fact that 

P + T ac Tp s u + = 6 a pu + (A. 6) 

which follows from Eq. (17), we obtain the following equation for <f>: 

= 2dJ3^. (A.7) 
Now we can rewrite the previous as follows: 

e / A((^e- / ) = <f>(e f Ae- f ) (A. 8) 
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The left hand side of the previous equation vanishes as a consequence of the 
self-duallity of the gauge field configuration. Hence, the Laplace operator 
acting on the product of cf> and vanishes except whenever vanishes. 
Taking for F = e~f the form (22), we conclude that (j> must be a linear 
combination of the following functions 

= £/ (A.9) 

(x - X( k) Y 

Therefore there are N linearly independent zero modes, in agreement with 
the index theorem . 

If we use the form Eq. (23), then again Eq. (A.9) gives solutions of 
the Dirac equation, but this time only N — 1 are linearly independent, as it 
should. 

In general given W = F -1 = and letting (for k = 1,...,N) be 
the position of its zeros, we can write the most general solution of the Dirac 
equation (Dtp^ = 0) as follows: 

= E d V°) Fl/2 d( ( x _ x (k)y ) u + (A.10) 

where d(k) are arbitrary complex constants. In obtaining Eq. (A.10) we have 
combined Eq. (11), Eq. (A. 4) and Eq. (A.9). These solutions coincide with 
those found by Grossman [24]. 

Concerning the normalization, one must again distinguish the cases given 
by formulas (22) and (23). In the first case, after some algebra one can show 
that choosing d^ = &ki is an orthonormal basis of the space of zero-modes 
for i = 1, . . . N. In the second case, one must form an unitary transformation 
of this basis, orthogonal to the vector d^ = pi which is trivial. Hence, there 
are only N-l elements of the basis. 

Now let us specialize to the dyon solution Eq.(24). The total topological 
charge diverges and there are infinitely many zeros of Fjjs- Hence, there is 
also an infinite number of zero-modes. A possible non-orthogonal basis of 
the space of solutions is then: 

u k (x) = ^( ^J^ K (A.11) 
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Since both F and W are periodic in time with period 2tt it is natural to 
consider the different independent zero modes to be independent Fourier 
components: 

i^\x + 2tt) = e i(32 *4>W(x ) (A.12) 

Formula (33) and (35) of Section 3 summarizes these results. Notice that 
f3 ranges from to 1, but the value of f3 = gives a trivial solution, a 
consequence of the same phenomenon which occurs for e~f of the form (23). 
All solutions are orthogonal to each other. 

To conclude let us write down the expression for the projector onto the 
zero-modes. Using a and a' for the 4-spinor indices and a, a' for the colour 
ones we may write 

= 2^,( — )&( — )) P ,frr„) oo ,( 7 07r^P)aa' (A-13) 

The equation follows from the expression of the zero-mode and by using Eq 
(12) for P+. 
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